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Abstract 

We give a simple physical derivation of the photon emission rate from the weakly coupled 
quark-gluon plasma connected with the collinear processes q — > 'yq and qq — > 7. The 
analysis is based on the light-cone path integral approach to the induced radiation. Our 
results agree with that by Arnold, Moore and Yaffe obtained using the real-time thermal 
perturbation theory. It is demonstrated that the solution of the AMY integral equation is 
nothing but the time-integrated Green's function of the light-cone path integral approach 
written in the momentum representation. 



1. In recent years much attention has been given to the problem of the photon radia- 
tion from the hot quark-gluon plasma (QGP). This interest is mostly motivated by the 
possibility of production of the QGP in the heavy-ion collisions that are under active 
investigation at RHIC, and will be studied soon at LHC. The thermal photon radiation is 
expected to be a good probe of formation of the QGP, since the photons can freely escape 
the plasma and give information about the plasma temperature, T. The photon emission 
rate per unit time and volume can be written as 

dN n B (u) 11/ t \ 

» = -J^ Imn S(-- k »- w 

where is the retarded photon polarization tensor, ub(uj) = (exp(u/T) — is the 
photon Bose-Einstein factor. Presently the calculation of the photon polarization tensor is 
usually performed within the thermal perturbative Quantum Chromodynamics (pQCD) 
methods. The application of the pQCD to description of the QGP includes the well-known 
hard thermal loop (HTL) resummation [1] which is based on separation of the hard quark 
and gluon modes with momentum p ~ T and soft collective modes with p ~ gT. The 
HTL approximation leads to the infrared finite leading order 0(g 2 ) contribution to the 
photon radiation rate connected with annihilation (qq — > 751) and Compton mechanism 
{qg^iq)) [2]. 

The analysis of the two-loop contribution (connected with the bremsstrahlung q — > 
7g and the induced annihilation qq — > 7 physical processes) to the polarization tensor 
demonstrated [3] that it turns out to be of the same order in the coupling constant 
as the leading order result. This fact is a consequence of the collinear nature of the 
bremsstrahlung and induced annihilation, which leads to the divergence of their cross 
sections in the massless limits. In the thermal bath this divergence is regularized by the 
thermal quark mass m* ft ~ gT, and as a result these processes give a contribution of the 
same order as the Born approximation. It turns out that the higher-order ladder diagrams 
give the contribution of the same order as the two- loop diagrams [4]. The calculation of 
the ladder diagrams has been performed by Arnold, Moore and Yaffe (AMY) [5]. The 
AMY [5] result can be written as 
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where q s is the charge of quark species type s, np{p) = l/(exp(p/T) + 1) is the thermal 
Fermi distribution, and f satisfies the integral equation 

i5Ef(p ± ,p, k) = 2 P± - g 2 C F T J ^C(q ± ) [f(p ± ,p, k) - f(p ± - q ± ,p, k)] . (3) 

Here, Cp = 4/3 is the quark Casimir, 5E = | A; (p>^_ + m 2 )/2(p + k)p\, m q is the asymptotic 
thermal quark quasiparticle mass, the kernel of the integral equation can be expressed 



2 



via the HTL resummed gluon propagator. It can be written in a simple form as C(q ± ) = 
ml/qKql + mjj) [6], where m D = gT[{N c + N F /2)/3]^ 2 is the Debye mass. The C(qJ 
includes the contribution of both the screened longitudinal gluon as well as the unscreened 
transverse gluon. 

Physically the ladder diagrams describe the Landau-Pomeranchuk-Migdal (LPM) sup- 
pression [7, 8] of the collinear bremsstrahlung/annihilation. In the last decade considerable 
progress has been made in understanding the LPM effect in QCD for fast partons with 
energy E ^> T propagating through the QGP. Two different approaches to the LPM ef- 
fect in this regime have been developed: the BDMPS formalism [9, 10] and the light-cone 
path integral (LCPI) formalism [11, 12, 13, 14] (for a review, see [15, 16]). In the LCPI 
approach the radiation rate is expressed via the Green's function of a two-dimensional 
Schrodinger equation with an imaginary potential. 

In the studies [9, 10, 11, 12] the QGP has been modeled by a static system of the Debye 
screened color centers. In the case of the photon radiation from the thermal partons one 
should treat on the same footing the partons from which the photon is radiated and 
the target partons. In the AMY analysis [5] it was done using the real time thermal 
perturbation theory. In the present paper we give an alternative simple physical derivation 
of the AMY result within the LCPI approach. We demonstrate that in the case of the 
weakly coupled QGP when one can separate the hard and soft modes the treatment of the 
LPM suppression for the photon emission from the thermal partons can be made similarly 
to the case of fast partons. We show that the solution of the AMY integral equation (3) 
is nothing but the time-integrated Green's function of the LCPI approach written in the 
momentum representation. 

Whether the model of the weakly coupled QGP, used in the present and AMY [5] 
analyses, is a reasonable approximation to the QGP which is produced in AA-collisions 
is still an open question. However, in principle, our derivation may have a broader range 
of applicability then the pQCD description of the QGP, if one calculates somehow the 
function C(q_|_) and the quark quasiparticle mass using the non-perturbative methods, or, 
say, extracts them fitting the experimental data on the jet quenching. 

2. We consider the radiation of the photons with uj > T which is controlled by the 
dynamics of the hard quark modes in the QGP. The average energy for hard partons 
~ 3T. In the weakly coupled QGP these hard partons undergo typically only small angle 
multiple scattering due to interaction with the random gluon fields at the momentum 
scale ~ gT <C T. The large angle scattering is a very rare process. The typical separation 
between the large angle scatterings is L/. a . ~ 1 / g A T. Using this scale, from the uncertainty 
relation AEAt ~ 1, one can obtain for the typical off-shellness of the hard modes 8m ~ 
(T/Lj.a.) 1 / 2 ~ g 2 T, which turns out to be much smaller than the thermal quasiparticle 
masses for quarks and gluons m*' 1 ~ gT. Note that the Li a . at the same time gives the 
scale at which the soft rescatterings generate the scattering angle about unity. It follows 
from the L-dependence of the mean squared momentum transfer q 2 ~ g 4 T 3 L connected 
with the random walk of hard partons due to small angle multiple scattering. 

The above estimates show that at the space scale considerably smaller than L La the 
hard parton trajectories are almost straight lines slightly distorted by soft rescatterings. 
Thus, at L <C Li a , the hard modes can be described as the quasiclassical plane waves 
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propagating in the soft random gluon field. It is important that in first approximation 
one can neglect the statistics effects in treating the small angle scattering. Indeed, the 
typical space scale for the soft gluon modes 1/gT is much larger than the typical separation 
between the hard partons ~ 1/T. For this reason from the point of view of the hard parton 
scattering the soft gluon field can viewed as a uniform field at the scale ~ 1/T. In the 
uniform field all hard partons of the same type scatter in the same way, and the small 
angle scattering leads simply to some shift of the distribution function in the momentum 
space. A well-known example of this phenomenon is shift of the Fermi-sphere for the 
strongly degenerated Fermi gas in a uniform electric field. Any statistics effects will be 
suppressed by some power of g. 

The above facts on the hard parton trajectories are important from the point of view of 
the treatment of the q — > and qq — > 7 processes in the pQGP. Let us first consider the 
photon bremsstrahlung. From the uncertainty relation one gets for the photon formation 
length Lf ~ 2(1 — x)E /xm 2 q 1 , where E is the initial quark energy, and x = uj/E is the 
photon fractional longitudinal momentum (we choose the z axis along the initial quark 
momentum). For u ~ E ~ T one obtains Lf ~ l/g 2 T. Thus, one sees that Lf <C L La . It 
means that at leading order in g the photon bremsstrahlung is not affected by the hard 
scattering. Since the statistics effects can be neglected for small angle scattering as well, 
the bremsstrahlung occurs in the same way as for hard partons in a soft random gluon 
field. The only effect of the thermal bath of hard partons is the trivial Pauli-blocking 
for the final quark, which, of course, must be taken into account. Due to the collinear 
character of the induced photon emission, in accounting for the Pauli-blocking one can 
take the final quark momentum collinear to that for the initial quark. In this case the 
contribution of the bremsstrahlung into the photon emission rate can be written as 

dN br d br r°° 2 dP' ' (p,k) 



dtdVdk k 2 (2nf 



E jf dp P 2 n F (p)[l ~ n F ( P ~ k)]6(p - k) q 2dL ' (4) 



where d br = 4iV c is the number of the quark and antiquark states, dP q q (p, k)/dkdL is 
the probability of the photon emission per unit length from a fast quark of type s in the 
random soft gluon field without the thermal bath. 

The formation length for the qq — > 7 mechanism is also much smaller than Li a _. It 
means that the photon is produced in annihilation of almost collinear quark-antiquark pair 
induced by the small angle rescatterings in the soft random gluon field. The contribution 
to the radiation rate can be easily obtained using the detailed balance principle. For the 
plasma size much smaller than the photon absorption length it gives 

dN an / 7 m— 1 dN a hs 



dtdVdk 1 WJ dtdVdk 

The photon absorption rate on the right-hand side of Eq. (5) should be evaluated for the 
thermolized photons. It can be written via the probability distribution per unit length for 
the 7 — > qq transition dP* ^(k, p) / dpdL, where p is the final quark momentum, similarly 



1 This estimate does not account for the LPM suppression. If the radiation rate is suppressed by a 
factor Slpm, the real in-medium formation length is Lj PM <~ SlpmL/. 
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to the formula (4). Then, using (5) one obtains 



dN„ 



dtdVdk (2tt) 



E 



dP s 

dpnF(p)np(k — p)6(k — p) 



(k,p) 



dpdL 



(6) 



where k — p is the antiquark momentum, d an = 2 is the number of the photon helicities. 

3. The probability distributions entering (4) and (6) for the case of the random soft 
gluon field can be evaluated in the LCPI approach [11] similarly to the case of the static 
model of the QGP which has been used in [11] only for simplicity. For this reason we 
limit ourselves to highlighting only the aspects connected with the generalization to the 
arbitrary random external gluon field. The starting point of the LCPI approach is the 
representation of the parton wave function in the form 



ip(t, r) = exp(-iE£)<i>(z, £, p) 



(7) 



where ^ = t — z, r = (z, p), and it is assumed that E » m q . To leading order in m q /E 



the z-dependence of the transverse wave function 
two-dimensional Schrodinger equation 



(z, £, p) for fixed £ is governed by the 



■ d(j)(z,p) ~ , 

% — — — = H(f)(z,p) , 



with the Hamiltonian 



H 



[( P± - gA ± ) 2 + m\ 
2fi q 



+ g(A -A z ). 



(8) 



(9) 



Here, the Schrodinger mass is p q = E, p ± is the operator of transverse momentum, 
= (A , — Aj_, — A z ) is the external gluon vector potential (A^ should be understood as 
TqAp, for clarity we omit the color generators). The argument £ of is omitted in (8), 
since the £ dependence emerges only via the boundary condition for <fi. For a random 
external potential without loss of generality we can set £ = 0. For this choice the vector 
potential in the Hamiltonian (9) must be taken at x 11 = (t,p,z). Thus the Hamiltonian 
(9) describes the evolution of the partonic wave function along the light-cone. Eqs. (7), 
(8) hold for each helicity state. The photon wave function can be written in a similar way 
(of course, in this case the interaction and mass terms do not appear). 

The matrix elements for the q — > '-fq and 7 — > qq transitions can be written in terms 
of the Green's function for the Hamiltonian (9) and its photon analogue. In the LCPI 
approach they are written in the Feynman path integral form. The quark Green's function 
reads 



K q (p 2 ,z 2 \p 1 ,z 1 



) = J £>pexp ji J 



dz 



2 1 



m q 
2/i, 



(10) 



where p = dp/dz, = (1, p, 1). Since |p| 1 one can replace x^ by = (1, 0j_, 1). As 
compared with Ref. [11] where the static potential with only A (r) nonzero component 
has been used, now it is replaced by A' = u^A^x) with x^ = (z, p, z). 
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Let us consider the q — > transition. The probability of the photon emission per 
unit length diagrammatically is given by the graph a of Fig. 1, where the longitudinal 
coordinate of one of the vertex is fixed. Note that due to the relation — f * = fq the quark 
Wilson lines in the path integral for quark propagators turns out to be replaced by the 
antiquark Wilson lines in the complex conjugate quark propagators. 

Similarly to the static potential one can perform in the diagram la averaging over 
the external potential at the level of integrands. Then the interaction with the external 
potential turns out to be transformed into the interaction between trajectories which 
depends only on the relative transverse distances between the trajectories. This fact 
allows one to perform some path integrations analytically, and transform the graph la 
to lb. (the details can be found in [14, 15, 16]), in which the only dynamical variable is 
the transverse separation between quark and photon p = p q — p 7 , and the antiquark is 
located at the center of mass of the g7-system. The propagator describing the evolution 
of a spurious gg7-system in graph b of Fig. 1 in the path integral form reads 



J Vpexp ji J 



dz 



M{x)p 2 



1 



<S>({p},z 1 ,z 2 ) , 



(11) 



where x = u/E is the photon fractional momentum, M(x) = Ex(l — x) is the reduced 
Schrodinger mass, Lf = 2E(1 —x)/xm 2 q is the photon formation length which comes from 
the mass terms in the quark propagators (10), and 

<f>({p},z u z 2 ) = ((expj-is J dz[A q (p q ,z) + A^(p q ,z)]'j)) (12) 

is the averaged Wilson line factor. The quark and antiquark transverse coordinate in 
(12) in therms of p are given by p q = xp, p q = 0. For nonzero p l 2 the propagator 
(11) is not gauge invariant without the Wilson lines connecting the quark trajectories 
at zi t 2- This fact, however, is not important since in the final formula for the radiation 
rate this propagator only appears for the point-like initial and final states, i.e., for the 
configurations shown in Fig. lb. The formula for the probability of the photon emission 
per unit length is similar to that for the static potential (we use here the fractional photon 
momentum x instead of k in (4)) 



dP s 
dxdL 



u 

2Re J dzg(x) [K(p x , 0|p 2 , z) - K^Pi, 0|p 2 , z )} 



(13) 



Pl,2- 



where K vac is the free propagator for $ = 1 2 . The factor g in (13) is the vertex operator 
which accumulates the summation over the particle helicities. It is given by 



iW = 5EnMA})v(,{A»(A-^ 



y / 2 



'_d_ 

dPx 



_d_ 

d Py 



(14) 



V(x,{\}) = 



-iq s 



2M(x) V 2x 



[A 7 (2 -x) + 2\ q x\ , 



(15) 



2 In the case of an infinite plasma the vacuum term in (13) gives zero contribution. However, it is 
necessary for selection of the induced radiation for a hard parton produced in the medium 
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where {A} is the set of helicities for the q — > jq transition, as in (2) q s is the quark charge. 
The vertex factor (14) corresponds to the averaged over the quark color states probability, 
the corresponding normalization factor 1/N C is canceled by the N c which comes from the 
summing over the colors in the quark loop in the diagram of Fig. lb. 

Taking the advantage of the fact that the correlation radius for soft modes is much 
smaller than the photon formation length, the Wilson line factor $ can be written in the 
exponential form 

(16) 



r f z2 

$({p},z 1 ,z 2 ) = exp -/ dzP(xp(z)) 

J zi 



where, to leading order in the coupling constant, P(p) diagrammatically is given by the 
graphs a, b, c of Fig. 2. The wavy lines in Fig. 2 correspond to the gluon correlator (the 
color indexes are omitted) 

G{x-y)=u»u v {(A»{x)A»{y))). (17) 

The parallel straight lines in Fig. 2 show the quark trajectories along the light-cone t — z, 
and the ^-integrations along trajectories are performed with a fixed position of one of the 
vertex. We neglect the contribution to the function P(p) from the higher order diagrams, 
some of the diagrams of the order g A are shown in Figs. 2d, e, f. This approximation 
assumes that the interaction of hard partons with soft modes should be sufficiently small 
at the longitudinal scale about the Debye screening radius. One can say, that neglecting 
the higher order contribution to P(p) corresponds to the dilute gas approximation for the 
gluon correlators. In this approximation one gets 

oo 

P(p) = g 2 C F J dz[G(z, ± z) - G(z, p, z)] . (18) 

— oo 

Here we have taken into account the effect of the quark/ ant iquark color generators, which 
leads to appearing of the quark Casimir factor for each gluon correlator. Note that the 
QED analogue of (18) has been introduced in [17] in the analysis of propagation of the 
relativistic positronium atoms in amorphous media. One can write P(p) as the momentum 
integral 

P(P) = j ^- 2 [l ~ exp(.pq ± )]P(q ± ) . (19) 

One can easily show that if one uses for the gluon correlators their HTL expressions the 
P(q ± ) can be expressed via the function C(q ± ) which enters (3) as 

P(qJ = g 2 C F TC(q ± ) . (20) 

Note that for the static model of the QGP the function P(p) reads 

P(p) = =^, (21) 
where n is the number density of the color centers, and u q q(p) is the dipole cross section 

*M = C T C F a s J d q± {ql + ml)2 ■ (22) 
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Here Ct is the color Casimir of the scattering center. It is interesting that numerically the 
P{p) evaluated in the static model turns out to be close to that calculated with the HTL 
formula for C(q ± ) in the region p < 1/mo which is important for the radiation processes 
in the QGP. 

For the interaction factor written in the form (16) the propagator (11) is, evidently, 
the Green's function for the Hamiltonian 

" = -2Fw(|) ! -^ ) + I7' < 23) 

The formulas for the radiation rate in terms of the solution of the Schrodinger equation 
with the Hamiltonian (23) with smooth boundary conditions which are convenient for 
numerical calculations can be found in [13, 18]. 

To establish the connection between the above formulas of the LCPI approach and 
those of AMY [5] we rewrite formula (13) in terms of the in-medium light-cone wave 
function of the spurious gg7-system \l/ m (x, p, {A}). For \ q > = \ q it reads 



y m (x,p,{\}) = V(x,{\}) (J^-t\ y ^j I dzlC(p,0\p',z) 



(24) 



p=o 



The analogue of (24) in the vacuum case after ^-integration gives the well-known light-cone 
wave function of the Fock component in the quark 



^(x, p, {A}) = V(x, {A}) expt-iX^rrigK^pmgx) , 
here K 1 is the Bessel function. In terms of \I/ m formula (13) reads 

dP ^-Re^V*(x,{\}) (A -a, A) * m (x,p,{\}) 



dxdL 



9p x 



| A I V-'/'.l' 

The function \I/ m satisfies the integral equation 



p=o 



with 



/f Mix) 
dzfC vac (p,0\p', z) = J exp 

— oo — oo 

In momentum representation it takes the form 



M(x)(p-p'Y \z 



2\z\ 



p_i_ -+- e 



2iM(x) f dq ± 



pl + e 2 J (2tt) 



(25) 



(26) 



* m (x, p, {A}) = V(x, p, {A}) - J dp'W(p, p')P(xp')* m (x, p', {A}) (27) 



(28) 



(29) 
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with e = xm q . Evidently the \l/ m can be written in the form 

* m (x, P± , {A}) = V(x, {\})[F X ( P± ) - tX,F y ( P± )} , (30) 
where now F is the solution to the integral equation 

F(P ± ) = ^ + ^ / ^P (qL )[F(p ± ) - F(p ± - sqj] . (31) 

Making use of (15) and (30) after summing over the helicities one can rewrite (26) in 
momentum representation as 

dP^ q _ 2q 2 a em [l + (1 - x) 2 } ^ r dp±_ p p ^ 



dxdL xM(x) J (2tt 

Similar calculations for 7 qq transition give 

dP% qq _ 2q 2 N c a em [x 2 + (1 - x) 2 }^ f dp 



dxdL M(x) Im i (27p P ^ F ' (33) 

where now x = E/u, M(x) = cux(l — x), and F satisfies the integral equation 

F <PJ = 5iTK + |t| / ^«U)P(PJ " F(P, - <U>] • (34) 

Substituting (32) and (33) into (4) and (6) after adding together the bremsstrahlung 
and annihilation contributions one obtains the total photon emission rate. It can be 
represented exactly in the form (2), (3) after the change of variables p ± — > rrp^ for 
the bremsstrahlung, if one uses the function f connected with F by the relations: f = 
— i4E(l — x)F for bremsstrahlung, and f = —i4ujx(l — x)F for annihilation contributions, 
respectively. The difference in the variables for bremsstrahlung is connected with the fact 
that we use the z-axis along the initial quark momentum (which is natural from the point 
of view of the light-cone wave function for q — > 75 transition), and in [5] it is along the 
photon momentum. In the collinear approximation both the choices are equivalent. 

4. The proposed derivation of the radiation rate is based on the picture of the weakly 
coupled plasma that is formally valid at g <^ 1. In reality, for the QGP produced in AA- 
collisions, when T < (2 — 3)T C , g ~ 1.5 — 2. It was proposed that for such temperatures 
the plasma is a strongly coupled QGP (sQGP) [19]. Nonetheless one can expect that even 
for the sQGP the intuitive physical picture of photon emission will be the same, and our 
formulas should be a good approximation (at least for sufficiently large photon energies) 
if one somehow calculates P{p) using nonperturbative methods, say, uses the lattice sim- 
ulation, or AdS/QCD. What is really crucial for our calculations is the approximation of 
small parton scattering angles at the scale of the photon formation length. To understand 
the situation with this approximation one can use the results of the analysis of the RHIC 
data on jet quenching which give the direct information about P(p), and allows one to 
estimate the p± broadening for the thermal partons. In terms of P(p) the L-dependence 
of the p± distribution can be written as [17] 

HP±) = J7^y J dpexp[ip ± p - P{p)L\ . (35) 
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For the purpose of the qualitative estimates one can use the model of the static Debye 
screened centers used in [18] for the analysis of the RHIC data on the nuclear modification 
factor Raa- Using the oscillator approximation of the dipole cross section a(p) ~ Cp 2 , one 
obtains (p 2 ± ) = 2CnL. Evaluation of C in the double gluon exchange approximation with 
a s frozen at the value a s as 0.7 at low Q 2 [20, 21] gives C ~ 0.5 for scattering on a quark. 
This approximation allows one to describe the RHIC data on the nuclear modification 
factor [18]. In the oscillator approximation the typical in- medium photon formation length 
which accounts for the LPM suppression reads Lj PM ~ — x) /Cnx 2 ] 1 ^ 2 [13]. Let us 
consider the final quark which has the smallest energy, at the scale ~ Lj PM one gets 
(P±)/P 2 z ~ 2[Cn/cuE 2 (l - xf} 1 ' 2 . For u ~ 6T and E ~ 2cu one obtains (p 2 ± )/p 2 z ~ 0.3. 
With increase of oo this ratio will be even smaller. In fact, the angles will be somewhat 
smaller if one includes the collisional energy losses which can give the contribution about 
20-30% for RHIC conditions. Thus, one sees that even though the typical angles are not 
very small, nonetheless, the collinear approximation seems to be not far from reality, at 
least for sufficiently large photon energies, which, in fact, are most interesting from the 
experimental point of view. In principle, the fact that the collinearity is not strongly 
violated for RHIC conditions does not contradict to the large value of g. It simply says 
that the qualitative estimates which just count the powers of g in the situation when the 
exact values of the coefficients are unknown, are not very reliable. 

5. In summary, we have proposed a simple alternative derivation of the photon emission 
rate in the pQGP using the LCPI approach [11]. The new derivation is interesting for 
several reasons. First, the LCPI formulation is convenient for the numerical calculations. 
One can use the methods previously developed for the case of fast partons [13, 18], which 
reduce the calculation of the radiation rate to solving the two-dimensional Schrodinger 
equation with smooth boundary conditions. Secondly, since the LCPI approach is formu- 
lated in the coordinate space, it can be helpful for clarifying the applicability limits of the 
formulas obtained for a uniform plasma when they are applied to the expanding finite-size 
plasma produced in AA-collisions. Say, the LCPI formulation can be used to estimate the 
role of the boundary effects connected with the initial stage of the QGP evolution when 
the radiation rate is the largest. Indeed, in this region, say, for the proper time r < 1 fm, 
probability of the photon emission in the partonic processes may be substantially reduced 
due to the restriction on \z\ < r in the integral on the right-hand side of (13). This is 
not possible in the AMY [5] formulation in the momentum representation. We leave the 
phenomenological applications of our results to further publications. 
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Figures 




a b 

Figure 1: The diagram representation of the x-spectrum for q — > ^q transition: (a) 
the starting form in terms of the quark and photon propagators, the left-directed lines 
correspond to complex conjugate propagators the vertical dotted line indicates the initial 
quark density matrix; (b) the final form corresponding to (13) in which the only dynamical 
variable is the transverse separation between photon and quark. 




a b c 




d e f 

Figure 2: The diagram representation of the function P(p) in (16): (a), (b), (c) the 
leading order contribution corresponding to (18); (d), (e), (f) some of the diagrams of the 
order of g A . 
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